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ditions  ai  frequencies  less  than  100  MHz  [4}.  This  is  because  mus¬ 
cle  is  primarily  resistive  below  100  MHz;  the  loss  tangent  (the  ratio 
of  the  loss  factor  to  the  dielectric  constant)  increases  from  2.1  ai 
100  MHz  to  5.8  at  13.56  MHz  for  dog  skeletal  muscle  [5].  Others 
have  used  phantoms  simulating  both  the  dielectric  constant  and  the 
conductivity  of  muscle  at  frequencies  below  100  MHz  (7]-|9),  and 
the  new  mixtures  described  may  also  be  considered  in  applications 
requiring  complete  simulation  of  the  dielectric  properties  of  mus¬ 
cle. 

We  do  not  report  thermal  properties  of  these  new  tissue-simu¬ 
lating  materials.  However,  other  phantom  materials  [30]  and  di¬ 
verse  biological  materials  [31]  with  high  water  content  have  heat 
capacities  and  thermal  conductivity  values  close  to  those  of  pure 
water,  and  we  expect  that  the  materials  described  here  will  as  well. 
For  precise  studies  involving  calorimetry,  the  heat  capacity  of  any 
phantom  should  be  measured  directly  by  the  investigator. 
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for  ■  microstrip  resonator  situated  on  an  isotropic  substrate  are  used 
to  validate  the  theory. 


I.  Introduction 

Microstrip  and  disk  resonators  are  among  key  important  ele¬ 
ments  in  the  design  of  microwave  and  millimeter-wave  frequency 
sources  and  filtering  networks  [1],  The  analysis  of  such  structures 
has  been  extensively  investigated  by  numerous  authors.  Different 
methods  based  on  various  approximations  [2]  to  the  quasi-static 
approach  13]-[4]  have  been  used  to  study  these  resonators.  Since 
these  analyses  are  based  on  approximations,  which  are  inaccurate 
for  high  frequency  applications,  quite  often  a  correction  factor  is 
needed  in  order  to  compensate  for  the  dispersion  effects.  More  ac¬ 
curate  theories  have  been  developed  and  implemented  to  study  the 
characteristics  of  the  microstrip  resonators,  as  documented  in  [5J— 
[6].  Up  until  now,  the  emphasis  of  the  resonator  studies  have  been 
directed  towards  structures  printed  on  isotropic  materials,  that  is, 
the  substrate  being  characterized  by  a  zero-rank  permittivity  ten¬ 
sor. 

In  this  paper,  a  full-wave  analysis  via  the  spectral  domain  ap¬ 
proach  is  used  to  study  the  characteristics  of  shielded  microstrip 
resonators  printed  on  anisotropic  substrates.  The  material  is  char¬ 
acterized  by  a  second  rank  permittivity  tensor  in  which  all  diagonal 
elements  can  be  distinct.  The  Green’s  function  of  the  structure  is 
derived  by  applying  the  resonance  conditions  and  by  solving  the 
fourth-order  differential  equations  in  the  Fourier  domain.  A  pro¬ 
cedure  based  on  Galerkin’s  method  is  then  used  to  form  the  char¬ 
acteristic  equation  from  which  the  resonant  frequency  of  the  reson¬ 
ator  is  numerically  obtained.  Data  for  a  resonator  printed  on  an 
isotropic  substrate,  which  is  treated  as  a  special  case  of  anisotropic, 
are  computed  and  then  compared  to  existing  data,  and  good  agree¬ 
ment  is  observed.  The  resonant  frequencies  of  the  microstrip  re¬ 
sonators  on  anisotropic  materials,  such  as  PTFE  cloth,  boron  ni¬ 
tride,  epsilam-10,  and  sapphire  are  then  computed. 


II.  Theory 

The  microstrip  resonator  structure  under  consideration  is  shown 
in  Fig.  1  along  with  the  coordinate  system  used.  A  rectangular  strip 
of  width  w  and  length  d  is  printed  on  one  side  of  the  anisotropic 
layer  situated  inside  a  shielded  housing.  For  simplicity  of  the  anal¬ 
ysis,  it  is  assumed  that  the  operating  frequency  is  below  the  cut¬ 
off  frequency  of  the  partially  loaded  waveguide.  The  substrate, 
which  extends  from  —  b/2  to  +b/ 2,  is  also  assumed  to  be  lossless 
with  thickness  h\.  The  distance  from  the  microstrip  resonator  to  the 
top  cover  is  h2.  In  general,  the  anisotropic  medium  is  characterized 
by  its  permittivity  tensor: 


Id  = 


0  0 
t„  0 
0  <„ 


and  |  mI  =  /UH 


(1) 


(b) 

Fig.  1.  Geometry  of  the  shielded  microslrip  resonator,  (a)  End  view,  (b) 
Top  view. 


The  Fourier  transform  of  any  quantity  is  defined  through  the  fol¬ 
lowing  relation: 

{*•  rb/2 

1  4(x,  y,  z)e'a7eilSl  dy  dz  (3) 

-•  J  -b/2 


with  a  and  B  being  the  Fourier  transform  variables.  When  (2a)  and 
(2b)  are  transformed  according  to  (3)  and  reduced,  a  set  of  coupled 
differential  equations  written  in  terms  of  the  tangential  electric  fields 
can  be  derived: 


d2£, 

+  y2£,  +  yt£,  =  0 
d*£.  r  „ 

+  u£y  =  0 


(4a) 

(4b) 


with  y2,  y4,  z2  and  z«  are  the  transformed  coefficients.  In  order  to 
find  the  solution  for  £}  or  £t,  (4a)  and  (4b)  are  then  decoupled 
yielding  a  fourth  order  differential  equation  for  cither  field.  Once 
the  tangential  field  solutions  are  found,  the  remaining  components 
in  the  anisotropic  region  are  then  derived  and  they  are  expressed 
in  their  generalized  forms  as 

£,  =  A*  sin  (y  +  x)  +  B *  sin  (y_x)  (5a) 

£y  =  Ey  A*  sin  (y.x)  +  E~  B*  sin  (y  x)  (5b) 

=  (7'wm.)'1  H* A*  cos  (y.x)  +  ( juipj  1  //,  B*  cos  (y.x) 

(5c) 


where  <„  and  are  the  frce-spacc  values. 

The  vector  wave  equations  for  the  components  of  the  electric  and 
magnetic  fields  within  the  anisotropic  layer,  which  cun  be  manip¬ 
ulated  from  Maxwell’s  equations,  arc  expressed  in  their  simplified 
forms  as 

(2a) 

(2b) 


ft,  -  (jwn„y'  H*  A*  cos  (7 ,  Jr)  +  '  ny  B*  cos  (7.x) 


(5d) 

with 

^  -  {(7 1 )?  .’..l/c, 

(6a) 

e:  .  rt'“| 

(6b) 

7*  (  <>.  1 

VX(VXE)  -  A’M  /•;  0 

VX(|«]  1  •  VXH)  kill  -  0 
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Hr  =  7±  ~  jfiE*  (6c) 

Ht  =  ~y  iEy  +j<*E?  (6d) 

where  -y±  can  be  obtained  by  solving  the  fourth  order  characteristic 
equation  through  the  method  described  in  [7]. 

The  fields  within  the  isotropic  region  2  can  be  derived  based  on 
the  superposition  of  TE  and  TM  fields.  The  transforms  of  the  scalar 
potentials,  which  are  used  to  find  the  field  components,  are  written 
below: 

d2S'k  . 

+  (yi>  V-  =  o.  (7) 

After  applying  the  boundary  conditions  at  the  top  cover  x  -  A, 
+  h2,  the  field  expressions  for  region  2  are, 

£,2  =  -  c-iB.)  sin  y2(hx  +  h2  -  x)  (8a) 

=  c35,  sin  y2(h,  +  h2  -  x)  (8b) 

ft, 2  =  {o*A„  +  CiB.)  cos  y2(h,  +  h2  -  x)  (8c) 

ft  a  ~  ciA*  4:05  7i(7»i  +  h2  -  x)  (8d) 

with  c,  to  c6  are  the  known  constants. 

At  the  air-anisotropic  layer  interface  jc  =  /i, ,  the  appropriate 
boundary  conditions  are  also  imposed  so  that  a  set  of  matrix  equa¬ 
tions  can  be  formed  yielding  the  expression  for  the  impedance 
Green’s  function: 


Ga(k„, «.  0) 

Gzy(k0.  a,  0) 

0) 

£t(a,  0) 

Gn(ka,  a,  0) 

G„(k0,  or,  0) 

7, (or,  0) 

£,(a,  0) 

with  the  matrix  elements  defined  as 

Ga  =  //,/det.  Ga  =  -//j/det 
Gy.  =  -//j/det,  Gyy  =  HJ det 


and  then,  substituting  7,  and  7,  into  the  impedance  Green's  function 
(9a),  and  after  taking  the  inner  products  with  7 *  and  7„  for  different 
indexes  i,  a  set  of  algebraic  equations  is  derived: 

M  N 

2  1)C„  +  2  2 )Dm  =  0  /  =  1.  2.  3,  •  •  •  N 

m»t  m - I 

M  N 

2  Tlm(2.  1)C„  +  2  7^(2.  2 )Dm  =  0  /  -  1 , 2.  3.  •  •  •  M 

*■  I  M-  I 

(11a) 

with 

n.0.  IK*.)  =  2  f  /„(or,  0)Gaik„  a ,  0)Jm(a,  0)  d0  (1  lb) 

m  -  I  J 
0 

m 

Tji  1,  2)(Jt.)  =  2  f  7*(a,  f3)<3w(*„,  a.  |3)7„( a.  13)  dp  (I  lc) 

I  J 
0 

09 

T„(2,  «(*.)  -  2  [  7„(a,  0)  Gy,(k0,  a.  0)JM(a.  0)  d0  (lid) 

H  ™  I  J 
0 

•  ’  09 

Tm{ 2.  2)(Jt0)  «  2  (  7„(a,  0)G„(k„,  a.  0 0)  d0.  (1  lc) 

H  “  l  J 
0 

The  right  hand  sides  of  (1  la)  have  been  eliminated  in  the  Galerkin 
process  through  the  application  of  Parscval’s  theorem.  The  simul¬ 
taneous  equations  are  then  solved  for  the  wave  number  k0  by  setting 
the  determinant  of  the  coefficient  matrix  equal  to  zero  and  search 
for  the  root  of  the  resulting  equation. 

For  the  dominant  mode,  J,,  and  J„  arc  chosen  to  have  the  fol¬ 
lowing  forms: 

/„(«,  13)  =  J,(a)J2(0)  (12a) 


det  =  HJh  -  H2H}  (9b) 

where 

ft,  =  UK2  +  {-{«**  +  ?,*'}/*„  (9c) 

//j  =  ?7*J  }/*, o  (9d) 

ft 3  =  +  {Ms**  -  Ms*- }/£,,  (9c) 

h.  =  *,*2  +  {-Ms**  +  M6*~'}/t,o  m 

*2  =  cot  (y2h2).  (9g) 

**  cot  (y./t,).  (9h) 

*  =  cot  (y.h,).  (9i) 


/„(«,  13)  =  73(a)7«(|3)  (12b) 

where  the  inverse  transforms  of  Jx,  J},  7Jt  and  74  are  given  by 

J,  =  sin  ((2ry/w)/[I.O  -  (2 y/w)2]l/2  (12c) 

J2  =  2z/d2  (1 2d) 

J,  =  1.0  -  (2y/w),r,/J  (12c) 

J,  ~  (2/d)  cos  (xz/d).  (12f) 


Note  that  the  form  of  Jx  and  Jy  is  identical  to  the  one  used  for 
computing  the  dispersive  characteristics  of  infinitely  long  micro¬ 
strip  lines 

III.  Results 


The  remaining  constants  in  (9c)-(9i),  £,  to  £n,  can  be  written  in 
terms  of  the  transformed  variables  as  well  as  the  medium  parame¬ 
ters. 

To  find  the  resonant  frequency  of  the  microstrip  resonator,  a  pro 
eedurc  based  on  the  Galerkin  technique  (6|  is  used,  by  first  ex¬ 
panding  the  unknown  currents  7,  and  7,  in  terms  of  known  basis 
functions 

M 

7,(r«,  0)  =  2  0)  ( 10a) 

m  “  I 

N 

Jiit.  0)  -  2  0)  (i()h) 


To  validate  the  theory,  the  resonant  frequency  of  a  shielded  mt- 
crostrip  resonator  printed  on  an  isotropic  substrate,  which  is  treated 
as  a  special  case  of  anisotropy,  is  computed  and  compared  to  pub¬ 
lished  data.  The  dimensions  of  the  structure  arc  b  =  155.0  mm,  h, 
12  7  mm,  h2  »  88.9  mm,  and  w  ==  20.0  mm.  The  substrate  is 
characterized  by  its  permittivity  tensor  as  <„  =  tr,  =  <„  «=  2  65 
Fig  2  shows  the  response  of  the  resonant  frequency  of  the  reso 
nator  with  respect  to  different  physical  lengths.  As  can  be  seen,  the 
resonator  may  resonate  anywhere  between  0.5  GHz  to  0.8  GHz 
range,  depending  on  the  chosen  dimensions.  Plotted  also  with  out 
data  are  the  results  reproduced  from  reference  [6],  and  an  excellent 
agreement  is  observed 
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Fig.  2.  Resonant  frequency  versus  the  resonator  length  with  =  t„  = 
=  2.65,  b  =  155.0  mm.  h,  =  12.7  mm,  h2  =  88.9  mm,  and  w  =  20.0 
mm. 


RESONANT  FREQUENCY  (OTl) 


Fig.  3.  Resonant  frequency  of  the  microstrip  resonator  printed  on  aniso¬ 
tropic  substrates  with  b  =  155  0  mm,  h,  =  12.7  mm,  h2  —  88.9  mm,  and 
w  =  20.0  mm. 

Next,  the  data  for  microstrip  resonators  printed  on  anisotropic 
substrates  are  computed.  For  the  biaxial  substrate,  namely  PTFE 
cloth,  the  medium  parameters  arc  c,,  =  2.45,  «ry  =  2.89,  and 
=  2.95.  For  the  uniaxial  cases,  they  arc  specified  as  =  3.4  and 
f>r  =  {:r  =  5.12,  f,,  =  10.3  and  <yv  =  t.,  =  13.0,  and  «„  =11.6 
and  <>r  =  =  9.4,  written  respectively  for  boron  nitride,  epsilam- 

10,  and  sapphire.  Fig.  3  shows  how  the  resonant  frequency  re¬ 
sponds  when  various  materials  arc  used.  Note  that  the  physical  pa¬ 
rameters  used  here  arc  the  same  as  for  those  previously  specified 
for  the  resonator  printed  on  the  isotropic  substrate.  The  results  show 
that  by  changing  the  medium  from  PTFE  cloth  to  boron  nitride, 
the  resonant  frequency  reduces  considerably,  and  this  is  in  contras! 
to  what  happens  when  the  medium  changes  from  epsilam-10  to  sap¬ 
phire.  In  the  second  case,  only  a  small  variation  in  resonant  fre¬ 
quency  is  observed  This  behavior  can  be  confirmed  by  using  trans¬ 
mission  line  theory,  where  the  resonant  frequencies  of  the 
resonators  arc  approximated  by  (X„/2 Since  the  effective 
dielectric  constant  of  an  infinitely  long  microstrip  line  on  cpsilam- 
10  is  very  close  to  that  computed  for  the  sapphire  substrate,  this 
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implies  that  the  resonant  frequencies  of  the  resonators  for  these  two 
media  are  just  about  the  same. 

IV.  Conclusion 

An  analysis  based  on  the  spectral  domain  method  applied  to  study 
shielded  microstrip  resonators  printed  on  anisotropic  materials  was 
presented.  The  anisotropic  layer  is  generally  specified  by  its  second 
rank  permittivity  tensor.  The  Green’s  function  for  the  structure  is 
obtained  through  a  fourth  order  D.  E.  formulation,  Galerkin’s 
method  testing  procedure  in  the  Fourier  domain  is  applied  to  form 
the  characteristic  equation  from  which  the  resonant  frequency  of 
the  resonator  is  numerically  obtained.  Data  on  the  resonant  fre¬ 
quency  of  resonators  primed  on  both  uniaxial  and  biaxial  substrates 
were  also  generated. 
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Analysis  of  a  Coupled  Slotlinc  on  a  Double-Layered 
Substrate  Containing  a  Magnetized  Ferrite 

Masahiro  Gcshiro  and  Tatsuo  Itoh 


Abstract—  Nonreciprocity  in  the  propagation  characteristics  of  the 
even  and  odd  modes  in  magnelized-fcrrite-loaded  double-layered  cou¬ 
pled  slotlines  is  studied.  The  analysis  is  based  on  Galcrkin’s  method 
applied  in  the  Fourier  transform  domain.  Numerical  results  are  pre¬ 
sented  for  various  values  of  structural  parameters.  As  a  result,  it  is 
found  that  the  waveguide  structures  studied  have  sufficient  nonreci¬ 
procity  in  propagation  constants  for  isolators  and  four  port  circula¬ 
tors. 

I.  Introduction 

Slotlincs  arc  well  suited  for  their  usage  in  nonreciprocal  ferrite 
devices  for  microwave  integrated  circuits  since  the  magnetic  field 
has  elliptical  polarization  (l|.  Recently,  several  analytical  methods 
for  propagation  characteristics  in  slotlines  ami  finlincs  containing 
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